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ABSTRACT 


The object of this report is to determine an underlying or basic 
theoretical statistical model for making probability inferences in regard 
to thunderstorm activity at Cape Kennedy, Florida. The negative binomial 
distribution with probability density function 
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is presented to represent the variation in thunderstorm events per day 
at Cape Kennedy. Statistical theory and methods are developed using the 
latest and most comprehensive thunderstorm data available. The conclu- 
sion is reached that the negative binomial distribution is the logical 
choice for an underlying model to represent thunderstorm events at Cape 
Kennedy, Florida. 


NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 




TECHNICAL MEMORANDUM X-53816 


A PROBABILITY DISTRIBUTION FOR THE NUMBER OF 
THUNDERSTORM EVENTS AT CAPE KENNEDY, FLORIDA 

SUMMARY 


This investigation was made to determine an underlying, or basic, 
theoretical statistical model for making probability inferences in 
regard to thunderstorm activity at Cape Kennedy, Florida. The negative 
binomial distribution with probability density function 


P(x) 


r(x + k) 
x! r(k) 


k x 

p (i - p) ; 


x = 0, 1 , 2, ... 

k > 0, 0 § p § 1 


is presented to represent the variation in thunderstorm events per day 
at Cape Kennedy. The statistical properties necessary for the applica- 
tion of the negative binomial are presented, and it is shown that these 
attributes are present in the distributions of thunderstorm events at 
Cape Kennedy. The latest and most comprehensive thunderstorm data avail- 
able are analyzed, and the conclusion is reached that the negative binomial 
distribution is the logical choice for an underlying model to represent 
thunderstorm events at Cape Kennedy, Florida. 


I. INTRODUCTION 


Statistical methods of analysis may be divided into two general 
categories, descriptive and analytical, both of which depend on the 
basic laws of probability. Descriptive methods reduce large amounts 
of data to a few meaningful "statistics" such as means and standard 
deviations. A theoretical statistical model (distribution function) is 
assumed for the observations, and analytical methods are used to deter- 
mine how well the empirical data fit this model; i. e . , the analytical 
procedures determine the "goodness of fit" between theory and observation. 

The purpose of this report is to determine an underlying, or basic, 
theoretical distribution for making probability inferences in regard to 
thunderstorm activity at Cape Kennedy, Florida. 



Thunderstorms are of primary concern in the design of launch vehicles, 
in the planning of space missions, and in launch operations at Cape 
Kennedy because of high winds, lightning hazard, and extreme turbulence 
associated with this atmospheric phenomenon. The combinations of environ- 
mental conditions, including unstable air with a relatively high moisture 
content, and some type of lifting action present during the summer months 
make Florida one of the major thunderstorm genesis areas over the entire 
earth. The negative binomial distribution is presented to represent the 
variation in thunderstorm events per day at Cape Kennedy. 

The author wishes to acknowledge the assistance of David Riggenbach, 
Terrestrial Environment Branch (R-AERO-YT) , for compiling the frequency 
distributions of thunderstorm events for Cape Kennedy from the National 
Weather Records Center data. 

II. STATISTICAL MODEL 


In practical statistics, a discrete probability law is required to 
describe events which seem to occur at random; for example, the arrivals 
of customers at a service point or the number of accidents and break- 
downs in a factory. It is common practice to assume that the frequencies 
of such events fit a Poisson distribution. However, the Poisson series 
requires the assumption that the probability of the event remains constant. 
This assumption implies that the variance of the distribution equals its 
mean. In reality, it is rarely true that the probability of the event 
remains constant. Any variation in the probability of the event, in 
particular, the tendency for one event to increase the probability of 
another, will increase the variance of the distribution in relation to 
the mean — which means a negative binomial distribution will better 
describe the data. A report by the weather observer of a thunderstorm 
is proof that the atmosphere is in a state of instability and conditions 
are present for the formation of further thunderstorm cells; i.e., the 
probability of the event is increasing. 

Let us consider the first application of the negative binomial 
probability distribution by Yule in 1910 [1], We will make an analogy 
between this application and the distribution of thunderstorms at Cape 
Kennedy. Suppose we have a population of people subjected to recurring 
exposures to a disease and that during an exposure each member of the 
population has an equal probability p of contracting the disease. After 
x exposures, the proportions who have contracted the disease 0, 1, 2, ... 
times will be given by 
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q , xpq 
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where q = 1 - p. The terms given by (1) are terms of the binomial series 
(q + p) x . If k unfavorable exposures to the disease are fatal to the 
individual, the proportion surviving after x exposures will be given by 
the first k terms of the binomial (q + p) x . The proportion dying during 
the xth exposure will be those who contracted the disease (k - 1) times 
in the first (x - 1) exposures and who contract it again during the xth 
exposure; i.e., it will be 


f x - l\ k-l x-k f x - l\ k x-k 

(js - ij* « • p = Vk - i)f « < 2 > 

and since deaths do not begin until the kth exposure, the proportion of 
deaths at the kth, (k + l)th ... exposure will be 


k 

P 



kq, k(k + 1) 



(3) 


which are successive terms in the expansion of p^(l - q)”^, a binomial 
with a negative index. Thus, the proportions of the original population 
dying during successive exposures are given by successive terms of the 
negative binomial distribution with the first deaths occurring at the 
kth exposure. 

Now, the probability of exactly x events (density function) is 
given by 


P(x) 


fx + k - l\ k x 
( k-l ) P * ‘ 


(4) 


Suppose in Yule's classic example we let the people exposed to the 
disease be analogous to the days in some month, say, June, being exposed 
to the synoptic condition favorable for the formation of thunderstorms at 
Cape Kennedy. Now, the number of deaths that result from exposure to the 
disease will be analogous to the number of thunderstorms that actually 
develop in June. Now, we have all the days in June subjected to recur- 
rent exposures of synoptic conditions favorable for the formation of 
thunderstorms. We must assume that each day in June that is exposed to 
the favorable synoptic conditions has an equal probability p of having 
a thunderstorm develop. This is a reasonable assumption. Continuing 
our analogy, the proportion of thunderstorms that develop at the kth, 

(k + l)th ... exposure will be given by (3), successive terms in the 
expansion of p^(l - q) - ^, a negative binomial whose density function is 
given by (4). 
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Thus, statistical theory indicates the negative binomial distribu- 
tion as the appropriate model for the distribution of thunderstorm 
activity at Cape Kennedy, Florida. 


III. ESTIMATION 


Numerous estimators for the parameters of the negative binomial 
distribution have been proposed. We have chosen to use the first two- 
moment method proposed by Cohen [2] . The negative binomial density 
function given by (4) may be written in terms of the gamma function as 


P(x) = 


r(x + k) 

x! r(k) 


k X 

p q , 


x=0, 1, 2, 


k > 0, 0 S p § 1. 
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The distribution function is given by 


F(x) 



T(x + k) 
x! r(k) 


k x 

p q , 


( 6 ) 


which gives the probability of obtaining a value of x less than or equal 
to some particular value of x, say x Q . 

Now, after some algebraic manipulation of Cohen's estimators, we 
have for the moment estimators of the parameters k and p 





P* = 


k + x 


(7) 


where x is the sample mean and s 2 is the 
The mean M of the negative binomial 


sample variance, 
distribution is given by 


M 





and the variance V is 
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The efficiency of estimating p and k by the method of moments is 
derived by Fisher [3], In terms of the parameters used here, the 
reciprocal of the efficiency is given by 


1 

E 


1 + 2 


(k + 2) 




2-3 


(k + 2) (k + 3) 


2-3-4 


(k + 2) (k 4- 3) (k + 4) 


IV. DATA SAMPLE 


According to standard United States weather observing procedure, a 
thunderstorm is reported whenever thunder is heard at the station. It 
is reported along with other atmospheric phenomena on the standard 
weather observer's form WBAN-10 when thunder is heard and ends 15 minutes 
after thunder is last heard. Notice that the standard definition of a 
thunderstorm may include multiple occurrences of thunderstorms. For 
this reason, we have chosen to use the term "thunderstorm event" as a 
more appropriate definition for our statistical analysis. 

Since reliable representative data concerning the number of thunder- 
storms actually passing over Cape Kennedy (or the launch site) is not 
available, it should be pointed out that the statistics presented in this 
paper are applicable only to an area surrounding Cape Kennedy defined by 
the distance at which thunder can be heard by the weather observer. An 
observer can hear thunder up to a radius of approximately 25 kilometers. 
The statistics are not appropriate for making probability inferences in 
regard to the number of thunderstorms that actually strike the launch 
pad , 


Also, the type of statistical analysis presented is useful primarily 
for the planning of missions rather than for application to operations. 
Statistics may be useful up to a few days before a mission. However, at 
this time the weather forecaster's predictions should be more accurate, 
and the transition is made from statistical inference to weather fore- 
casting dependent upon the synoptic situation prevailing a few days 
before the mission. 

The data sample used was produced by ESSA, National Weather Records 
Center, Asheville, North Carolina, under government order number H-76789 
for the Terrestrial Environment Branch, Aerospace Environment Division, 
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and is the latest and most comprehensive thunderstorm data available for 
Cape Kennedy, Florida. The period of record is January 1957 through 
December 1967. 


V. ANALYSIS 


The negative binomial and Poisson distributions were tried as 
prospective models for thunderstorm events at Cape Kennedy. Table 1 
summarizes observed frequencies of days that experienced x thunderstorm 
events for all months, and for the spring, summer, and fall seasons at 
Cape Kennedy. Table la gives the relative frequency of occurrence of 
days that experienced at least one thunderstorm event at Cape Kennedy 
for the same reference periods. 

Theoretical summaries of the months and seasons that experience 
significant thunderstorm activity at Cape Kennedy are given in tables 2 
through 12. In all cases, the sample variance was significantly greater 
than the sample mean, indicating the negative binomial distribution as 
the appropriate model. The Kolmogorov-Smirnov statistic was used for 
a ’’goodness of fit” test. 

Notations used in tables 1 through 12 are as follows; 

x = the number of thunderstorm events per day 

f = the observed number of days during the 11 -year period 

of record that experienced x thunderstorm events 

r.f. = the relative frequency of occurrence of x thunder- 
storm events 

F = the observed distribution function 

o 

f e = the expected frequencies using the negative binomial 

distribution 

F(x) = the negative binomial distribution function 

x = the sample mean 

s 2 = the sample variance 

k* , p* = parameter estimators of the negative binomial 
distribution 

n = sample size 
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= the tabulated Kolmogorov-Smirnov statistic for 
sample size n and rejection level a = 0.05 

|F - F(x) | = the maximum absolute difference between the 

observed distribution function and the negative 
binomial distribution function. 

Conditional probabilities are also included in the tables. 

Consider the month of June (table 5) as an example. There were 
40 days out of 330 days (11 years of Junes) that had exactly two thunder- 
storm events. This gives a relative frequency (probability) of occur- 
rence of 0.121 of having exactly two thunderstorm events during any day 
in June. The observed distribution function (F Q ) gives a probability of 
0.921 of having two or less thunderstorm events during any day in June, or 
a probability of (1 - 0.921 = 0.079) of having more than two thunder- 
storm events during any day in June. The negative binomial distribution 
predicts 36.9 days in June that will experience exactly two thunderstorm 
events and the probability (F(x)) is 0.928 of having two or less thunder- 
storm events during any day in June, or a probability of (1 - 0.928 = 
0.072) of having more than two thunderstorm events during any day in 
June. The agreement between theory and observation is very good. Com- 
paring F 0 with F(x) shows a maximum absolute difference in the distribu- 
tion functions of 0.017 occurring at x = 0. Since the Kolmogorov-Smirnov 
statistic 


D 


n 

a 


= D 


330 
O. 05 


is equal to 0.075, the value of 0.017 is not sufficiently large to 
reject the hypothesis at the 5 percent rejection level that this sample 
can be fitted by a negative binomial distribution. 

T r, e conditional probabilities are computed from the theoretical 
frequencies (f e ) by using a double summation technique due to 0. E. 
Smith.* 

The tabulation on page 8 is an example of this technique using 
the month of June (see table 5). Each element in the second summation 


/V 

Chief, Terrestrial Environment Branch, Aerospace Environment Division, 
Aero-Astrodynamics Laboratory, Marshall Space Flight Center, Alabama. 
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(ZZ) is divided by the appropriate top element in each column as indicated 
in order to obtain the conditional probabilities; i.e., in each column 
under conditional probabilities, given i thunderstorm events (i = 1, 2, 3, 
the probability of having k additional thunderstorm events (k = 0, 

1 , 2, . . . ) is given by 


P(i + k | i) = 


(i 3- k)th element 
(i)th element 


( 8 ) 


For example, for i = 2; given two thunderstorm events on any day in June, 
what is the probability of having two additional thunderstorm events 
(k = 2) on that day in June? From equation (8), 


p ,, | _ 4th element _ 12.5 

^ ' ' 2nd element 96.2 


0.130. 


Also, given four thunderstorm events on any day in June (i = 4), the 
probability of having one additional thunderstorm event (k = 1) on that 
same day in June is 0.304. 


VI. CONCLUSIONS 


There are many advantages in the use of a theoretical statistical 
model for predicting a variable such as thunderstorm events at Cape 
Kennedy. Once sufficient representative samples have been collected and 
analyzed and the validity of the theory is established by an appropriate 
statistical test, the theoretical model becomes "deterministic" and may 
be applied universally to the variable under consideration. Another 
advantage of theory over empirical statistics is the use of the accept- 
able theoretical function for making probability inferences concerning 
values of the variable outside of the range of observation. It is often 
desired to make predictions relating to these "never observed" values, 
and the theoretical approach permits one to do so. It should be pointed 
out that no theoretical function can explain all observations for which 
it is the proposed model. Some areas of non-agreement must occur between 
theory and observation. These areas should be considered as expected 
deviations of the observations from the "fitted" theoretical curve. 

The physical properties necessary for the application of the nega- 
tive binomial distribution have been shown to be present in our experi- 
ment concerning the number of thunderstorm events at Cape Kennedy. In 
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all the samples considered, the sample variance exceeded the sample mean, 
indicating the negative binomial as the appropriate model. Our comparison 
with Yule's classic application of the negative binomial substantiates its 
validity to represent the number of thunderstorm events at Cape Kennedy, 

In all 11 samples considered, the conclusion was reached that the 
negative binomial distribution could not be rejected at the 5 percent 
rejection level as being an acceptable model for the number of thunder- 
storm events per day at Cape Kennedy. Furthermore, the negative binomial 
gave a "better" fit in all cases than the Poisson distribution. 

Using statistical theory and methods, we have demonstrated that the 
negative binomial distribution is the logical choice for an underlying 
model to represent thunderstorm events at Cape Kennedy, Florida. 
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TABLE 1. Frequencies of the Observed Number of Days that Experienced x Thunderstorm Events 
at Gape Kennedy, Florida for the 11 -year Period of Record January 1957 through 
December 1967 * 
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TABLE 2. March-Negative Binomial Distribution for Thunderstorm 
Events at Cape Kennedy, Florida 


X 

f 

o 

r.f . 

F n 

f 

F(x) 

Conditional Probability 







i 







0 

308 

.902 

.902 

305.4 

. 896 

1 2 3 4 

l 

20 

.059 

.961 

25.5 

.970 

1 

2 

9 

.026 

.987 

6.7 

.990 

.277 1 

3 

3 

.009 

.996 

2.2 

.996 

.078 .281 1 

4 

1 

.003 

1.000 

.8 

.999 

.016 .059 .211 1 

X = 

0.150 

s 2 - 

0.268 

k* = 0,189 

p* = 0.558 n = 341 


Goodness of 

fit: 

D n = 0.074 

a 







o 

i 

X 

'w' 

1 “ 0 

.009 

TABLE 3. 

April- 

Negative Binomial Distribution for Thunderstorm 
Events at Cape Kennedy, Florida 

X 

f 

r.f. 

F 

f 

F(x) 

Conditional Probability 


o 



e 


i 







0 

299 

.906 

.906 

295.9 

.897 

1 2 3 

i 

18 

.055 

.961 

25.3 

.973 

1 

2 

10 

.030 

.991 

6,1 

.992 

.226 1 

3 

3 

.009 

1.000 

1.8 

.997 

.042 .186 1 

X = 

0.142 

s 2 = 

0.237 

k* - 0.214 

p* = 0.600 n = 330 


Goodness of 

fit; 

D a = °- 075 







|F 0 - F(x) 

1 = o 

.013. 
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TABLE 4: May-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 


X 

f 

o 

r.f . 

F 

o 

f 

6 

F(x) 


Conditional 

Probability 







i 

0 

266 

.779 

.780 

262.6 

.770 

i 

2 3 

4 5 6 

l 

43 

.126 

.906 

52.4 

.924 

i 



2 

25 

.073 

.979 

16.6 

.972 

.339 

1 


3 

3 

.009 

.988 

5.9 

.989 

.120 

. 354 1 


4 

3 

.009 

.997 

2.2 

.996 

.041 

.122 .345 

1 

5 

0 

.000 

.997 

.9 

.998 

.013 

.037 .106 

.306 1 

6 

1 

.003 

1.000 

.3 

.999 

.003 

.007 .021 

.061 .200 1 


X = 

= 0.352 

s 2 = 0. 

621 k* 

= 0.460 p* 

= 0.567 n = 

341 



Goodness 

of fit: 

4= o 

.074 





F q - F(x) | = 0.017 


TABLE 5: June-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 



f 

r f 

T? 

f 

F(x) 


Conditional 

Probability 

X 

o 

X . i. , 

r 

o 

e 

i 

0 

187 

.567 

.567 

181.5 

.550 

i 

2 3 

4 5 6 

i 

77 

.233 

.800 

87.7 

.816 

i 



2 

40 

.121 

.921 

36.9 

.928 

.394 

1 


3 

17 

.052 

.973 

14.7 

.972 

,147 

.373 1 


4 

6 

.018 

.991 

5.7 

.989 

.051 

.130 .348 

1 

5 

2 

. 006 

.997 

2.2 

.996 

.016 

.040 .106 

.304 1 

6 

1 

.003 

1.000 

.8 

.999 

.003 

.008 .022 

.064 .211 1 


X = 

= 0.752 

s 2 = 1. 

169 k* 

= 1.354 p* 

= 0.643 n = 

= 330 



Goodness 

of fit: 

D n = 0, 

a 

.075 





F q - F(x) | = 0.017. 
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TABLE 6: July-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 


X 

f 

o 

r.f. 

F 

o 

f e 

F(x) 

Conditional Probability 

i 

0 

177 

.519 

.519 

166.2 

.487 

i 

2 

3 4 5 

l 

80 

.234 

.753 

99.4 

.779 

i 



2 

47 

.138 

.891 

45.4 

.912 

.399 

1 


3 

26 

.076 

.967 

18.6 

.967 

.143 

.357 

1 

4 

9 

.026 

.993 

7.2 

.988 

.044 

,110 

.307 1 

5 

2 

.006 

1. 000 

2.7 

.996 

.009 

,023 

.066 .214 1 


x = 0.874 

s 2 = 1.277 

k* = 

1.893 

P' f = 

0.684 

n = 341 


Goodness of fit; D n = 0,074 

a 

|F q - F (x) | = 0.032 


TABLE 7: August-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 


X 

f 

o 

r.f. 

F 

o 

f 

e 

F(x) 

Conditional 

Probability 

i i 1 

0 

185 

.542 

.542 

180.2 

.528 

1 

2 

3 4 5 

i 

89 

.261 

.803 

92.2 

.799 

1 



2 

30 

.088 

.891 

40.5 

.918 

.399 

1 


3 

24 

.070 

.961 

16.9 

.967 

.146 

.366, 

1 

4 

10 

.029 

.990 

6.8 

.987 

.046 

.116 

.316 1 

' 5 

3 

.009 

1.000 

2.7 

.995 

! .010 

.026 

.070 .221 1 


x = 0. 

,809 

s 2 = 1.280 

JU 

k = 

1.391 

* 

P - 

0.632 

n = 341 


Goodness of fit: D n = 0.074 

a 

|f q - F(x)| = 0.026 
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TABLE 8: September-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 


X 

f 

o 

r.f. 

F 

0 

f 

e 

F(x) 

Cond itional Pr ob ab i 1 i ty 

i 

0 

228 

.691 

.691 

219.2 

.664 

1 2 

3 4 

l 

54 

.164 

.855 

73.1 

.886 

1 


2 

33 

. 100 

.955 

24.8 

.961 

.316 1 


3 

12 

.036 

.991 

8.5 

.987 

.089 .283 

1 

4 

3 

.009 

1 . 000 

.2.9 

.995 

.018 .057 

.203 1 

X 

= 0.509 

s 2 

= 0.777 

k* = 

0.967 

p* = 0.655 

n = 330 


Goodness of fit: D n = 0.075 

a 

|F q - F (x) | = 0.031 


TABLE 9: October-Negative Binomial Distribution for Thunderstorm 

Events at Cape Kennedy, Florida 


X 

f 

o 

r.f. 

F 

0 

f 

e 

F(x) 

Conditional Probability 

i 

0 

311 

.911 

.911 

307.7 

.902 

12 3 

l 

17 

.050 

.961 

24.2 

.973 

1 

2 

9 

. 026 

.987 

6.1 

.991 

. 235 1 

3 

4 

.012 

1.000 

1.9 

.997 

.045 .192 1 

X = 

0.138 

s 2 

= 0.242 

k* = 

0.182 

p* = 0.570 n = 341 


Goodness of fit; D^ = 0.074 

I F - F(x) I = 0.011 
i o ' 1 
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TABLE 10: Spring (March, April, May) -Negative Binomial 

Distribution for Thunderstorm Events at Cape Kennedy, Florida 


X 

f o 

r.f. 

F 

o 

f 

e 

F(x) 

Conditional Probability 


0 

873 

.863 

.863 

863.6 

.853 

■ 

2 

3 4 

5 6 

i 

81 

.080 

.943 

103.7 

.956 

i 




2 

44 

.043 

.986 

29.2 

.985 

.312 

1 



3 

9 

.009 

.995 

9.8 

.994 

. 106 

. 339 

1 


4 

4 

.004 

.999 

3.5 

,998 

.035 

.113 

. 335 1 


5 

0 

.000 

.999 

1.3 

. 999 

.011 

. 034 

.101 .303 

1 

6 

1 

.001 

1.000 

.5 

1 . 000 

.002 

.007 

.022 .066 

.217 1 


x = 0.215 s 2 = 0.386 k* = 0.271 p* = 0.557 n = 1012 


Goodness of fit; d” = 0.043 

|F - F (x) 1 = 0.013 


TABLE 11; Summer (June, July, August) -Negative Binomial 
Distribution for Thunderstorm Events at Cape Kennedy, Florida 


X 

f 

o 

r.f. 

F 

o 

f 

e 

F(x) 

Conditional Probability 

i 

0 

549 

.542 

.542 

527.8 

.522 

i 

2 

3 

4 5 6 

l 

246 

.243 

.785 

279.6 

.798 

i 




2 

117 

.116 

.901 

122.7 

.919 

.404 

1 



3 

67 

.066 

.967 

50.1 

.969 

.153 

.379 

1 


4 

25 

.025 

.992 

19.7 

.988 

.054 

.133 

.352 

1 

5 

7 

.007 

.999 

7.6 

.995 

.017 

.041 

.108 

.307 1 

6 

1 

.001 

1.000 

2.9 

.998 

.004 

.009 

.023 

.067 .216 1 


x = 0. 

812 

s 2 = 1. 

245 k* 

= 1.523 p' ! 

' = 0. 

652 

n = 1012 


Goodness 

of fit: 

D n = 0. 

a 

043 






F 0 - F (x) | = 0.021 
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TABLE 12: Fall (September, October, November) -Negative Binomial 

Distribution For Thunderstorm Events at Cape Kennedy, Florida 


X 

f 

o 

r.f . 

Fo 

f 

e 

F(x) 

Conditional Probability 

i 

0 

860 

.859 

.859 

845.2 

. 844 

12 3 4 

l 

77 

.077 

.936 

109.5 

.954 

1 

2 

45 

.045 

.981 

30.6 

.984 

.286 1 

3 

16 

.016 

.997 

10.1 

.994 

.080 .281 1 

4 

3 

.003 

1.000 

3. 6 

.998 

,017 .058 .208 1 


x = 0. 

227 

s 2 = 0.397 

k* = 

0.302 

P* = 0.571 n = 1001 


Goodness of fit: = 0.043 

|F - F(x) | = 0.018 
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